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SUMMARY 

ABSTRACT.  The  US  Army  Concepts  Analysis  Agency  (CAA)  is  responsible  for  modeling  the 
employment  of  tactical  nuclear  weapons  at  the  theater  level.  Currently  available  simulations  that 
model  the  exchange  of  tactical  nuclear  weapons  require  much  time  to  set  up  and  run.  CAA  has 
developed  a  simple  probability  model  based  on  alternating  renewal  processes  to  provide  a  quick 
estimate  of  the  probability  that  a  relocatable  combat  unit  can  be  acquired  and  engaged  by  nuclear 
weapons,  without  a  loss  in  accuracy  compared  to  the  detailed  simulation.  The  model  may  be  applied 
to  either  tactical  or  strategic  relocatable  targets. 

THE  RESEARCH  SPONSOR  was  the  Director.  US  Army  Concepts  Analysis  Agency. 

THE  OBJECTIVE  OF  THE  RESEARCH  was  to  develop  an  analytic  model  for  the  acquisition  and 
movement  of  relocatable  nuclear  targets,  in  order  to  estimate  the  probabilities  that  target  units  may 
be  acquired,  be  available  for  fire  planning,  and  be  available  for  engagement  with  a  nuclear  weapon. 

THE  MAIN  ASSUMPTIONS  used  in  this  research  were: 

(1)  For  any  target  unit,  the  time  to  acquire  it,  the  time  that  it  is  retained,  the  time  that  it  stays 
in  place  and  the  time  that  it  is  moving  are  mutually  independent  random  variables  with  stationary 
probability  distributions. 

(2)  The  acquisition  process  and  movement  process  for  any  given  unit  are  independent. 

(3)  Units  are  acquired  and  move  independently  of  each  other  during  the  short  timeframe  of 
interest  (generally  less  than  12  hours). 

THE  BASIC  APPROACH  used  in  this  research  was  to  model  the  acquisition  process  and  the 
movement  process  for  each  unit  as  independent,  stationary  alternating  renewal  processes. 

THE  PRINCIPAL  FINDING  of  the  research  is  that  is  is  possible  to  develop  a  model  of  the 
acquisition  and  movement  state  of  relocatable  nuclear  targets  using  an  alternating  renewal  process 
representation. 

THE  RESEARCH  WAS  PERFORMED  BY  MAJ  Mark  A.  Youngren. 

COMMENTS  AND  QUESTIONS  may  be  sent  to  the  Director,  US  Army  Concepts  Analysis  Agency. 
ATTN:  CSCA-RQN,  8120  Woodmont  Avenue.  Bethesda,  MD  20814-2797. 
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Modeling  Ike  Acquisition  and  Engagement  of  Relocatable  Nuclear  Targets 


Introduction 

The  US  Army  Concepts  Analysis  Agency  (CAA)  performs  analyses  for  the  Department  of  the 
Army  relating  to  engagement  of  armed  forces  at  the  theater  level,  including  the  tactical  employment 
of  nuclear  weapons.  When  we  examine  issues  relating  to  tactical  nuclear  warfare,  we  need  Lo  be  able 
to  determine  the  probability  that  we  can  acquire  potential  nuclear  targets,  retain  them  long  enough 
to  plan  the  nuclear  fires,  and  have  the  target  still  be  in  place  at  the  time  that  detonation  occurs.  In 
order  to  alter  the  course  of  the  battle  to  meet  theater-level  or  strategic  objectives,  nuclear  weapons 
use  may  be  constrained  to  achieve  a  specific  purpose  within  an  appropriate  period  of  time.  For 
example,  weapons  may  not  be  fired  at  potential  targets  as  they  are  acquired;  there  may  be 
significant  delays  between  the  time  the  fires  are  planned  and  the  time  of  detonation. 

CAA  currently  uses  the  Nuclear  Fire  Planning  and  Assessment  Model  III  ( NUFAM  III) 
(Schuetze  and  Albrecht  [1980]),  a  detailed,  stochastic  simulation  model  of  a  two-sided  tactical 
nuclear  exchange,  to  support  analyses  of  theater  nuclear  issues.  The  model  is  very  useful  in  providing 
a  representation  of  the  effects  of  a  postulated  nuclear  exchange,  but  it  is  very  time-  and  resource- 
intensive  to  set  up,  run,  and  analyze  the  output.  This  situation  generated  a  desire  for  a  simple 
analytic  model  that  can  be  used  to  quickly  estimate  the  probability  of  engagement  of  potential 
nuclear  targets  without  having  to  set  up  a  series  of  simulation  runs. 

Our  solution  to  this  problem  is  a  renewal-based  probability  model  of  acquisition  and 
engagement  of  relocatable  nuclear  targets  at  the  tactical  or  strategic  level  that  can  be  used  to 
quickly  estimate  the  probability  that  targets  can  be  acquired,  planned  for  fire,  and  be  available  to  be 
hit  by  a  nuclear  weapon.  NUFAM  III  currently  explicitly  simulates  independent  movement  and 
acquisition  times  using  the  same  movement  and  acquisition  distributions  throughout  the  simulation 
-  thus.  NUFAM  is  generating  realizations  of  independent,  stationary  alternating  renewal  processes, 
and  the  probability  model  described  herein  provides  a  direct  alternative  to  the  estimation  process 
internal  to  NUFAM.  Of  course.  NUFAM  III  also  models  processes  not  addressed  in  this  paper,  such 
as  weapon  allocation  and  damage  assessment.  However,  use  of  this  renewal  model  in  NUFAM 
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[II  will  greatly  increase  the  efficiency  of  the  simulation,  and  use  outside  of  NUFAM  will  permit  a 
quick,  identical  estimate  of  the  probabilities  of  acquisition  and  engagement  without  having  to 
exercise  the  NUFAM  simulation. 

Modeling  Theater- Level  Tactical  Nuclear  Warfare 

Modeling  nuclear  weapons  at  a  theater  level  poses  several  problems,  brought  about  by  an 
asymmetry  in  the  resolution  of  our  representation  of  theater  combat  and  the  resolution  required  for 
modeling  tactical  nuclear  effects.  Combat  at  a  theater  level  involves  large  areas  of  space  and  vast 
numbers  of  combat  units,  support  forces,  and  militarily  significant  entities.  Limitations  brought 
about  by  available  computer  resources  (even  using  supercomputers)  and.  more  importantly,  our 
capability  to  understand  and  properly  represent  such  combat  in  a  computer  model,  require 
significant  simplification.  The  most  important  simplifications  in  theater  level  modeling  involve  the 
relatively  low  level  of  resolution  of  the  models,  and  the  use  of  deterministic  (as  opposed  to 
stochastic)  processes  and  data.  We  rely  upon  these  low-resolution  combat  models  or  scenarios  to 
provide  the  initial  situation  (force  arrays,  strengths,  etc.)  to  our  models  of  tactical  nuclear  warfare. 

Several  models  of  conventional  warfare  exist  at  the  theater  level.  The  model  used  at  CAA  is 
called  the  Force  Evaluation  Model  (FORCEM).  Like  most  theater-level  models  and  scenarios. 
FORCEM  simplifies  by  representing  combat  forces  at  the  division  and  higher  level  and  representing 
time  in  12-hour  time  steps.  Our  models  of  tactical  nuclear  warfare  such  as  NUFAM  III.  on  the  other 
hand,  require  a  much  higher  degree  of  resolution  in  time  and  space.  Combat  units  may  be  targeted 
at  the  company  level  and  artillery  units  at  the  level  of  individual  missile  launchers.  These  units  may 
move  and  be  acquired  multiple  times  within  12  hours.  The  period  during  which  nuclear  weapons  are 
planned  and  fired  is  modeled  on  a  time  scale  of  minutes  to  hours,  well  within  a  12-hour  time  period. 
As  a  result,  we  are  uncertain  about  the  precise  location,  movement  and  acquisition  of  these  higher 
resolution  units,  and  we  rely  on  probability  models  to  describe  our  uncertainty  regarding  them. 
Specifically,  we  need  to  represent  the  effects  of  acquisition  and  movement  of  potential  target  units 
within  a  period  of  time  (e.g.,  12  hours)  during  which  the  opposing  force  structures  remain 
approximately  the  same  (as  represented  in  theater-level  combat  models). 

Ideally,  we  would  like  to  represent  the  detailed  acquisition  and  movement  processes, 
interdependent  with  conventional  combat,  the  environment,  etc.  in  detailed  simulation  or  analytic 
probability  models.  As  we  have  explained  previously,  this  level  of  detail  is  not  practical  at  the 
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theater  level.  What  we  do.  instead,  is  to  run  high  resolution  models  of  small  force  sizes  and  areas  to 
derive  estimates  of  the  distributions  of  the  acquisition  and  movement  times  of  interest.  These 
distributions  are  averaged  over  stochastic  environments,  terrain,  etc.  to  derive  summary  estimated 
distributions  for  the  acquisition  and  movement  times  of  interest.  These  distributions  will  hold  for  at 
least  the  12-hour  period  within  the  conventional  model  or  scenario  which  defines  the  scenario,  force 
structure,  sensor  capabilities,  etc.  for  these  high  resolution  mddels.  In  practice,  these  averaged 
distributions  are  used  to  describe  our  uncertainty  about  acquisition  and  movement  over  longer 
periods  of  time  (24  to  36  hours).  Models  and  scenarios  other  than  FORCEM  may  have  different 
resolution  in  time  but  the  same  principles  apply. 

The  issue  of  primary  interest  in  this  paper  is  how  we  use  this  summary  information  to 
adequately  represent  the  acquisition  and  movement  processes  during  the  time  when  tactical  nuclear 
weapons  are  planned  and  fired.  NUFAM  HI  uses  these  distributions  to  generate  an  explicit  sequence 
of  acquisition  and  movement  events  for  every  unit  represented  in  the  model.  We  will  show  how  this 
explicit  simulation  of  acquisition  and  movement  can  be  eliminated  without  a  loss  in  relative 
modeling  fidelity. 

The  Target  Acquisition  Process 

Detecting  target  units  with  sufficient  accuracy  to  plan  for  nuclear  fires  is  the  process  of  target 
acquisition.  Targets  are  combat  units  or  major  items  of  equipment  such  as  mobile  missile  launchers 
that  are  considered  for  engagement  by  nuclear  weapons.  Relocatable  targets  are  targets  which  have 
the  capability  to  move  during  the  scenario  of  interest  (although  they  may  or  may  not  retain  mission 
capability  during  movement).  As  a  consequence  of  this  movement  capability,  targets  do  not  remain 
acquired  indefinitely  (unless  they  can  be  tracked  indefinitely  once  acquired);  at  some  time,  they 
move,  and  the  acquisition  is  no  longer  valid.  Even  if  a  tracking  capability  exists,  there  is  a 
probability  that  such  tracking  will  be  lost  over  time. 

The  outcome  of  the  target  acquisition  process  is  a  changing  acquisition  list.  A  target  unit  is 
acquired  when  it  is  detected  by  a  sensor,  identified  as  a  target,  and  placed  on  the  acquisition  list.  A 
target  unit  may  be  dropped  from  the  list  either  due  to  a  negative  sensor  report,  ( i . e . .  we  no  longer 
detect  its  presence),  or  it  may  be  dropped  after  some  period  of  time  when  the  acquisition 
information  cannot  be  updated.  Any  given  target  unit  will  alternate  between  two  states:  acquired 
(retained  on  the  list)  or  not  acquired.  A  target  acquisition  process  is  therefore  a  temporal  series  of 
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such  acquisition  states.  The  time  to  acquisition.  Ta ,  is  the  time  it  takes  to  acquire  a  target  once  any 
previous  acquisition  has  been  dropped:  the  tune  of  retention.  Tr,  is  the  time  a  target  is  retained  on 
t  he  acquisition  list. 

Targets  that  are  to  be  engaged  using  conventional  weapons  are  generally  fired  upon  soon  after 
acquisition.  Nuclear  targets  differ  from  conventional  as  they  are  planned  for  specific  purposes 
dictated  by  the  overall  tactical  and/or  strategic  situation.  As  a  result,  they  are  not  normally  engaged 
as  they  are  acquired:  rather,  nuclear  fires  are  directed  at  targets  that  are  acquired  and  perceived  to 
be  in  place  at  the  time  the  weapons  are  approved  for  fire. 

Figure  1  illustrates  a  representative  acquisition  sequence  for  a  relocatable  target  unit.  Once  the 
unit  has  been  dropped  from  the  list,  it  is  immediately  subject  to  being  reacquired.  We  expect  that 
nuclear  weapons  use  will  occur  after  the  conventional  battle  has  been  underway  for  some  time:  thus, 
we  are  interested  in  the  acquisition  probabifties  at  some  point  in  time  well  after  the  acquisition 
process  has  begun.  We  assume  that  the  times  to  acquisition  { Ta  }  are  independent  and  identically 
distributed  (iid)  and  make  the  same  assumption  about  the  times  of  retention  { Tr  } .  For  relocatable 
targets,  we  can  approximate  the  target  acquisition  process  as  an  alternating  renewal  process  of 
indefinite  length.  Both  the  time  that  the  target  is  dropped  from  the  list  and  the  time  that  the  target 
is  acquired  are  renewal  points  of  this  alternating  renewal  process. 

Retained  Tr 

Not  acquired  Ta 

Time  - ^ - 

Acquisition 

Figure  1.  Possible  Target  Acquisition  Sequence 

The  Fire  Planning  and  Engagement  Sequence 

Enemy  units  are  identified  as  potential  nuclear  targets  as  they  are  acquired.  If  the  scenario  calls 
for  the  use  of  tactical  nuclear  weapons,  the  model  must  be  able  to  represent  employment  against 
specific  types  of  targets  that  may  be  found  during  a  limited  period  of  time.  Because  of  the 
constraints  associated  with  nuclear  weapons  employment,  each  nuclear  fire  mission  must  be  carefully- 
planned  before  orders  are  given  to  fire  the  weapons.  The  fire  planning  process  will  normally 
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terminate  at  the  beginning  of  some  specified  time  interval  of  length  L  within  which  weapons  may  be 
fired,  which  we  refer  to  as  the  fire  execution  period.  We  identify  the  time  that  the  fire  planning 
process  ends  as  tP.  A  o.rget  unit  that  may  be  planned  will  have  had  a  recent  acquisition  time  A 
prior  to  tP.  The  nuclear  round  will  detonate  at  some  time  Tj  during  the  fire  execution  cycle: 
generally,  this  time  is  specified  during  the  planning  but  is  not  known  a  prion  when  modeling  the 
engagement  sequence  (Figure  2). 


j  Fire  planning 

Fire  execution  period 

) 

1 

I - ^ 

Time 

A 

tp 

Figure  2.  Acquisition.  Planning,  and  Detonation  Events  for  an  Acquired  Target 

The  important  thing  to  note  is  that  tP  is  chosen  independently  of  the  target  acquisition  aim 
movement  sequences  of  any  potential  target  unit.  As  a  result,  the  time  tp  may  be  viewed  as  a 
random  entry  point  into  these  processes.  Following  the  terminology  common  in  reliability,  we  call 
the  time  from  the  last  transition  to  our  random  entry  point  the  age  and  the  time  from  our  random 
entry  point  to  the  next  transition  the  residual  life.  An  example  using  the  acquisition  sequence  is 
shown  in  Figure  3. 


Retained  Age 

Residual  life 

Not  Acquired 

Time 

tP 


figure  3.  Example  Snowing  Age  and  Residual  Life  of  the  Retention  Time  Distribution 
The  Movement  Stale  of  the  Relocatable  Target 

The  target  unit  may  be  in  one  of  two  alternating  states  with  respect  to  movement:  it  can  be  in 
the  move  state  (moving),  or  it  cm  be  in  the  stay  state  (stationary).  We  define  the  random  variable  S 
to  represent  the  length  of  time  that  a  target  is  stationary  and  the  random  variable  M  to  indicate  the 
length  of  time  that  it  is  moving.  We  assume  that  all  S  and  M  are  mutually  independent  and 
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represent  the  unit  movement  as  an  alternating  renewal  process.  If  we  examine  the  subunit  state  at 
any  arbitrary  point  in  time,  it  has  a  probability  of  being  in  the  stay  state  (stationary)  equal  to  ps.ay- 

_  E[SJ  m 

p»<«»  E[S]  +  E[M]  ’  1  J 

where  the  notation  E[X]  denotes  the  expectation  of  the  random  variable  X. 

A  target  can  be  acquired  and  observed  while  either  stationary  or  moving,  so  Ta  and  Tr  may  be 
dependent  upon  S  and  M.  However,  modeling  acquisition  at  the  theater  level  requires  aggregation  of 
acquisition  data  representing  the  net  times  to  acquisition  of  more  than  10.000  potential  relocatable 
target  units  using  all  of  the  sensors  deployed  by  the  opposing  force.  The  only  theater-level  data 
available  at  the  present  time  is  generated  by  a  separate  target  acquisition  model.  TADER  (Penn  and 
Bauman  [1987]),  run  at  CAA.  Given  many  characteristics  of  the  sensors,  target  units,  and  battlefield 
scenario,  to  include  the  expectations  of  the  stay  and  move  times  S  and  M.  TADER  output  for  each 
unit  is  used  to  estimate  the  parameters  of  the  distributions  of  Ta  and  Tr.  The  effect  of  unit 
movement  is  averaged  out  in  the  TADER  calculations  and  is  not  recoverable.  As  a  result,  the  only 
distributional  data  available  to  support  the  determination  of  the  joint  distribution  of  Ta,  Tr.  S.  and 
M  are  the  univariate  distributions.  Therefore,  in  our  model,  we  assume  that  if  the  fixed  parameters 
of  the  distributions  of  Ta.  Tr,  S.  and  M  are  known,  then  the  variables  Ta.  Tr,  S.  and  M  are 
mutually  independent.  This  assumption  of  independence,  as  is  usually  the  case,  is  a  simplistic 
approximation  to  complex,  poorly  understood  dependencies. 

The  assumption  of  independence  between  the  processes  of  acquisition  and  movement  is  explicitly 
or  implicitly  present  in  all  models  and  studies  examining  nuclear  exchanges  at  the  theater  level  (to 
include  Nl'FAM  III)  of  which  the  author  is  aware.  At  the  high  level  of  aggregation  present  in 
theater  analyses,  even  more  simplistic  models  are  frequently  used.  For  example,  the  issue  of  target 
retention  is  frequently  ignored,  and  the  movement  of  units  is  either  ignored  or  represented  in  a  less 
realistic  manner.  It  is  our  intent  to  provide  a  simple  alternative  to  detailed  simulations,  such  as 
Nl’FAM  III.  that  p-  rapiu  calculation  of  the  probabilities  of  acquisition  and  engagement  of 

relocatable  nuclea.  -  ■  <’  j.  without  loss  of  precision.  If  the  model  presented  herein  is  used  in  lieu  of 
a  model  making  s;milar  -  more  simplified  assumptions,  this  purpose  is  achieved. 

Representing  Ike  Target  Acquisition  Process 

When  constructing  a  model  of  the  acquisition  and  engagement  of  relocatable  nuclear  targets,  it 
is  not  necessary  to  explicitly  represent  the  target  acquisition  and  movement  processes  in  a  detailed 
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simulation.  Since  these  processes  can  be  represented  as  alternating  renewal  processes,  a  well-known 
result  of  renewal  theory  provides  us  the  probability  of  acquisition,  pa<;?.  at  any  random  point  of 
time: 


p  ilrj 

Pae’  E[  Tfl  ]  +  E[  Tr  ]  ■ 


[-] 


The  time  to  acquisition.  Ta,  is  normally  derived  from  detailed  sensor  simulation  models.  The 
TADER  model  provides  for  each  unit  a  probability  pa  that  it  was  acquired  by  any  part  of  the  sensor 
systems  on  a  side  during  a  short  interval  of  time  6,.  These  values  pa  are  fit  to  a  standard  "glimpse" 
type  model  of  target  acquisition,  where  there  is  a  probability  pa  that  the  unit  is  detected  (acquired) 
during  any  time  interval  <5,.  Thus  the  probability  of  detecting  the  unit  on  the  Ath  such  glimpse  is 
pa(  1  —  Pa)*  —  \  A—  1.2 . This  is  geometric  with: 


w  * 

P[  Ta  >  t  ]  =  (  1— pa  ,  t  >  0;  E[Ta]  =  £, 


[3] 


where  [xj  denotes  the  largest  integer  less  than  or  equal  to  x. 


The  length  of  time  that  the  target  unit  is  retained  on  our  acquired  target  list.  Tr.  will  be 
dependent  upon  our  assumptions  about  the  intelligence  collection  process  and  the  units  represented 
in  the  scenario.  When  the  target  unit  cannot  be  observed  after  acquisition,  we  use  a  constant 
retention  time  equal  to  the  expectation  of  the  residual  stay  time:  that  is.  Tr  =  E[  Y,  j.  When  the 
target  unit  is  subject  to  periodic  observations  subsequent  to  the  acquisition  time,  the  period  between 
observations  is  assumed  to  be  the  same  as  for  the  acquisitions  (  6t  ). 


'he  probability  of  contini"  ;  to  detect  the  target  unit  during  a  single  glimpse,  given  that  if 
was  detected  on  the  previous  gi  pse,  is  pr.  then  the  probability  that  the  target  acquisition  is  lost 
during  a  single  glimpse  is  (  1  —  pr  ).  and  the  retention  time  Tr  =  d>t  Kr.  where  Kr  is  distributed  as 
geometric(  1  —  pr  ).  Obviously,  pr  will  be  considerably  higher  than  the  probability  of  detection 
given  that  the  unit  /as  not  detected  previously.  We  handle  the  continuous  observation  cases  by 
using  the  same  glimpse  model  with  a  very  short  interval  6t.  Methods  used  at  CAA  for  estimating  Tr 
are  discussed  below. 


The  probability  that  the  unit  detected  at  time  t  will  remain  on  the  target  list  until  some  time  s 
is  simply  the  probability  that  the  time  remaining  until  the  next  transition  point  (dropping  the  unit 
from  the  list)  after  (  will  be  greater  than  (  s  —  t  ).  If  we  enter  a  renewal  process  at  some  random 
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point,  we  will  be  more  likely  to  enter  at  a  point  where  there  is  a  long  time  between  transitions 
rather  than  an  average  length  of  time.  Thus  the  distribution  of  the  residual  life  (time  remaining 
until  the  next  state  transition)  is  not  the  same  as  the  underlying  interarrival  distribution. 
Fortunately,  the  distribution  for  the  residual  life  has  been  worked  out  (cf.  Ross  [1983],  pg.  68).  If  we 
let  Yr  denote  the  residual  life  of  the  retention  time,  renewal  theory  tells  us  that 


P[  Yr  >  y]  =  1 


1 

E[Tr] 


P[  Tr  >  u  ]  du  . 


[4] 


Combining  this  probability  with  pac?  given  above,  the  probability  that  a  unit  is  acquired  at  time  l 
and  remains  on  the  list  for  s  additional  hours  is 


P[  on  the  acquisition  list  at  time  i  and  retained  from  t  to  t  -f-  s  ] 

E[  Tr  ]  -  J’=0  FTr(u)  du 

E[  Ta  ]  +  E[  Tr  ]  ’  [5] 

where  F T  (u)  denotes  P[  Tr  >  u  ].  For  additional  details  on  implementing  alternating  renewal 
A  r 

processes,  see  Youngren  [1989]. 

The  parameters  of  the  distributions  of  Ta.  Tr.  S.  and  M  are  determined  based  on  the  theater- 
level  scenario  during  the  time  period  (in  the  case  of  FORCEM.  a  12-hour  time  period)  that  includes 
the  time  tP.  We  would  expect  that  these  parameters  will  change  over  time:  however,  these  changes 
are  generally  gradual  in  our  scenarios,  and  we  are  only  interested  in  the  process  for  a  short  period 
(less  than  12  hours)  around  the  time  tp.  As  a  result,  we  assume  that  these  parameters  remain  fixed. 

An  Example 

Suppose  that  the  probability  pa  that  a  target  unit  is  detected  during  an  interval  <5,  is  equal  to 
0.2  with  fi,  =  1  hour.  Then  E[  T0  ]  =  =  5  hours. 

Suppose  that  the  stay  time  S  is  distributed  as  uniform  (  10  min.,  40  min.  )  and  the  move  time 
M  is  distributed  as  uniform  (  90  min.,  120  min.  ).  Then  E[  S  ]  =  25  min.  and  E[  M  ]  =  105  min. 
From  this,  we  can  calculate  the  probability  that  the  unit  is  stationary  at  any  random  point  in  time. 

Pjtatj' 

~  25  +°105  =  °'192' 


8 


CAA-RP-89-6 


We  can  also  calculate  the  distribution  of  the  residual  stay  time,  Yj: 


P(  \ >  y  )  -  '  E|sj 


u  —  0 


P[  S  >  u  ]  du 


0  <  y  <  10  min. 

10  min.  <  y  <  40  min. 

C40 

Tiie  expectation  of  the  residual  stay  time,  E[  Y3  ]  =  P[  Y,  >  v  ]  dv 

J  j i=o 


=  f5(  SO  -  4y  +  fj  ) 


-  14  min. 


Determining  the  Probability  that  a  Unit  is  Available  for  Fire 


In  order  for  a  unit  to  be  available  for  fire  by  nuclear  weapons,  it  must,  be  acquired  and  retained 
as  a  target  at  least  until  the  scheduled  firing  time.  The  target  is  evaluated  at  some  known  planning 
time  tp;  during  the  planning  period,  the  detonation  will  be  scheduled  to  occur  at  some  time  Td 
within  the  fire  execution  period.  Thus  a  target  that  is  on  the  acquisition  list  at  time  tp  must  be 
retained  at  least  until  Td;  this  occurs  when  the  residual  time  on  the  acquisition  list.  Yr.  is  greater 
than  or  equal  to  (  Td  —  tp  )  (Figure  2).  We  define  the  probability  P[unit  is  acquired  and  retained 
until  detonation]  as  pavai/.  Thus  pavall  —  P[  Yr  >  Td  —  tp  |  acquired  at  tp  ]  •  P[  acquired  at  tp  ]. 


We  calculate  p avall  using  the  residual  life  of  the  retention  time.  In  order  to  make  the 
development  complete,  we  define  R  as  the  maximum  permitted  value  for  Tr  and  L  as  the  maximum 
permitted  value  for  (Td  — tp);  0  <  R  <  oo;  0  <  L  <  30.  For  convenience  we  use  the  notation 
•‘P[  X  =  x  ]”  to  denote  the  probability  density  function  f^(x)  of  a  continuous  random  variable. 
Thus  (conditioned  on  being  on  the  acquisition  list  at  time  tp)  we  have: 


P[\r>Td—  tp|Td  — tp_  t]  —  L  —  ,  Fj  (u) 

M  -*■  r]  u=0  1  r 


du 


=  0 


0  <  t  <  R. 
t  >  R, 


and 


P[  Yr  >  Tj-tp  ]  =  P[  Yr  >  Td— tp  I  Td  — tp  =  t  ]  P[  Td-tp=  t  ]  dt. 


t=o 


Thus 


P«*«i<=  Yr  >  Td  —  tp  |  acquired  at  time  tp  ]  -  P[  acquired  at  tp 


IL  {  1  -  EIT71  L.  du  }  ' 


E[Tr 


E[  Ta  ]  +  E[  Tr  ] 

where  D  =  min  {  R,  L  };  0  <  Td  — tp  <  L  <  dc:  0  <  Yr  <  R  <  oo;  Tr  >  0;  tp  >  0. 


[6] 
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In  some  cases  (such  as  a  and  b  below),  we  also  require  that  the  unit  be  stationary  at  the  last 
time  that  it  was  observed.  When  this  occurs.  P  ^  is  multiplied  by  the  probability  that  the  unit 
was  stationary,  p,tay- 

Determining  the  Probability  of  Hitting  a  Relocatable  Nuclear  Target 

A  unit  can  be  hit  if  it  is  available  for  fire  and  is  stationary  at  time  Tj  at  the  place  where  it  was 
last  observed.  We  do  not  consider  the  availability  of  a  suitable  weapon,  in  range,  that  can  engage 
the  target  unit,  nor  do  we  assess  the  effects  of  the  detonation,  taking  into  account  delivery  system 
accuracy,  etc.  The  calculation  of  the  probability  of  hit  is  designed  to  tell  us  what  our  opportunity  for 
hitting  a  relocatable  target  might  be.  The  precise  calculations  of  the  probability  of  hit  depend  upon 
the  capabilities  of  the  sensor  system(s)  to  observe  the  target  after  acquisition.  Four  cases  that  have 
been  used  at  CAA  are  discussed  below. 

Case  a.  No  Capability  Exists  to  Observe  the  Target  after  Acquisition 

This  case  represents  the  situation  when  there  is  no  opportunity  for  any  sensor  to  observe  the 
target  between  the  acquisition  time  and  the  end  of  fire  planning.  This  may  occur  when  a  potential 
target  is  observed  by  some  sort  of  sensor  (e.g.,  aerial  reconnaissance)  with  a  single  mission  over  the 
area,  when  the  period  between  observation  opportunities  is  long  with  respect  to  the  planning  and 
firing  process,  or  when  the  target  is  able  to  conceal  itself  to  the  extent  that  subsequent  sensor 
missions  cannot  verify  either  the  presence  or  absence  of  the  target.  In  this  case,  the  target  remains 
on  the  acquisition  list  until  there  is  a  perception  that  the  information  is  too  "old'’  to  plan  a  nuclear 
fire.  The  time  interval  that  defines  ’‘old”  will  be  a  matter  of  judgment  and/or  doctrine,  and  is 
generally  based  on  an  estimate  of  the  time  that  the  acquired  unit  will  remain  in  place:  e.g.,  the 
expected  residual  life  of  the  stay  time. 

In  order  for  the  acquisition  to  be  useful,  the  target  unit  must  have  been  stationary  at  the  time  of 
acquisition.  The  joint  probability  that  the  unit  is  on  the  acquisition  list  at  the  planning  time  tP  and 
was  stationary  at  the  acquisition  time  A  is: 

P[  unit  on  the  acquisition  list  at  tp  and  stationary  at  A  ] 

—  D  .  D  _  1 Jr]  .  _ E[SJ _ 

Pstay  E[  Ta  ]  +  E[  Tr  ]  E[  S  ]  +  E[  M  ]  ’  li 
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The  probability  that  the  unit  is  available  for  fire  is: 


p4!.a!i  =  P[  unit  retained  until  Td,  on  the  acquisition  list  at  tp,  and  stationary  at  A 
=  P(  Yr  >  Td— tp  |  on  acquisition  list  at  tp  and  stationary  at  A  ] 

•  P[  unit  on  the  acquisition  list  at  tP  ]  •  P[  unit  stationary  at  A  ]  .  thus 


fD 


J  t=o 


P[  Yr  >  Td— tp  |  Td  —  tp  —  t  ]  P[  Td  —  tp  —  t  ]  dt 


E[T, 


E[  S 


E[  Ta  ]  +  E[  Tr  J  E(  S  ]  +  E[  M  ]  • 


where  D  =  min  {  R.  L  };  0  <  Td  — tP  <  L  <  cc:  0  <  Yr  <  R  <  oc:  Tr  >  0:  A  >  0. 


[8] 


Finally,  the  probability  of  hit.  pA  is  the  joint  probability  that  the  unit  is  available  for  fire  and 
will  still  be  stationary  at  time  Td  at  the  location  at  which  it  was  acquired.  The  acquisition  time  A  is 
not  known,  but  we  know  that  the  difference  (  t P  —  A  )  is  distributed  as  the  age  of  the  distribution 
of  the  time  of  acquisition,  Yr  (the  age  and  the  residual  life  are  identically  distributed),  independent 
of  the  distribution  of  Td—tp.  Because  of  independence,  pAit  is  calculated  as: 


=  P[  Yr  >  Td  —  tp  |  on  acquisition  list  at  tp  ]  •  P[  unit  on  the  acquisition  list  at  tp  ] 
•  P[  Ys  >  Td  —  A  |  stationary  at  A  ]  •  P[  unit  stationary  at  A  ]  . 


=  r  r  pi  Yr  >  t,  -  tp  i  Td — tP = t  j  p[  Td— tP = t  ] 

L  j  t=o 


dt 


E[  Tr] 


E[  Ta  ]  +  E[  Tr  ] 

I  [  P[Yj>(Td~tp)  +  (tp  — A)  |  Td— tp  =  t  and  (tp  — A)  =  y]  P[Td  — tp  =  t]  P[tp  — A]=y]  dt  dy 
J  v=o J t=o 

E[S] 


E[  S  ]  +  E[  M  ]  * 


[9] 


where  D*  =  min{  L,  Su  }:  S«  =  max{  S  };  0  <  S  <  S «  <  oc. 


Example  Case  a 

For  this  example,  we  set  the  retention  time  Tr  ■=  E[  Y»  ]  =  14  minutes,  a  constant;  thus 
P[Tr>t]  =  1  for  0  <  t  <  E[  Yj  ];  0  for  t  >  E[  Y,  j.  Obviously  E[  Tr  ]  =  E[  Y,  ]  =  14  minutes. 
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Recalling  that  E(  Ti7  ]  =  5  hours, 


14 


=  0.045.  The  residual  retention  time  Yr  is  distributed  as: 


(5-60)  4-  14 

P[  Yr  >  }'  ]  =  1  —  wir-j  [  P[  Tr  >  u  ]  du  =  1  —  0  <  y  <  R;  0  for  y  >  R. 

rj  J  u  —  o  R 

where  R  =  max{  Tr  }  =  14  min. 


We  assume  that  the  detonation  time  is  equally  likely  to  occur  during  the  fire  execution  period 
from  tP  to  tp-f-L,  with  L=60  min.  Thus  (  Td  —  tp  )  ~  uniform  (  0,  60  min.  ) 


Pa  fail " 


r* 


t=0 


P[  Yr  >  Td-tp  |  Td-tp  =  t  ]  P[  Td  —  tP  =  t  ]  dt 


E[  T, 


E[S 


E[  Ta  }  +  E[  Tr  ]  E[  S  ]  +  Ef  M  ] 


1=0 


r  i _ l 

1  1  R 


11  .  J, 

R J  60 


dt 


•  Pac,  •  PsUy  -  0.117  ■  0.045  •  0.192  =  0.001 


To  calculate  p engage,  we  recall  that  (  tp  —  A  )  is  distributed  as  the  age  of  the  retention  time 
(identical  to  the  distribution  of  the  residual  life  Yr).  Thus 


Vtngagt  —  Pa„al;  ‘ 
°  I'D 


V=0  J 


t=0 


P[Y»>(Td—  tp)  +  (tp  —  A)  |  Td  — tp  =  t  and  (tp  — A)  =  y]  P[Trf  — tp  =  t]  P(tp-A  =  y]  dt  dv. 


After  some  messy  algebra, 

P engage  =  0.001  •  0.138  ~  0.0001  . 

Case  b.  The  Target  is  Observed  Periodically  after  Acquisition 

This  case  represents  the  situation  when  there  is  an  opportunity  for  at  least  one  sensor  to  observe 
the  target  periodically  (but  not  continuously)  after  acquisition.  This  may  occur  when  a  potential 
target  is  observed  by  some  sort  of  sensor  which  can  observe  the  area  at  periodic  intervals.  In  this 
case,  we  assume  that  the  target  remains  on  the  acquisition  list  until  a  negative  sensor  report  is 
received.  Note  that  this  may  occur  either  because  a  unit  has  moved  or  the  sensor  failed  to  detect  the 
unit  on  a  subsequent  observation.  The  probability  of  detecting  the  target  unit  on  a  relook  will 
normally  be  considerably  higher  than  the  initial  probability  of  detection. 
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We  define  a  time  T,  as  the  last  time  that  the  target  unit  was  observed  prior  to  tp;  thus 
T,  <  t p  <  Td.  The  analysis  from  this  point  onward  is  straightforward;  except  for  the  distribution  of 
(tp  —  T(),  all  of  the  equations  in  the  previous  section  hold  with  T,  replacing  A. 


P  tvatl 


=rr 

J  t=0 


P[  Yr  >  td  — tp  |  Td  —  tp  —  t  J  P[  Td  — tp  —  t  ]  dt 


E[  Tr  ] 


_E[S] 


E[  Ta  1  +  E[  Tr  ]  E[  S  ]  +  E[  M  ]  ' 


[10] 


Finally,  the  probability  of  hit,  pAl(,  is  the  joint  probability  that  the  unit  was  stationary  at  the 
time  of  acquisition  A,  it  is  on  the  acquisition  list  at  time  tp,  it  will  be  retained  as  a  target  at  least 
until  time  Td,  and  it  will  still  be  stationary  at  time  Td  at  the  location  at  which  it  was  last  observed. 
The  time  T,  is  not  known,  but  we  know  that  time  tp  occurs  at  a  random  point  during  the  interval 
between  observations.  We  can  regard  the  interval  {tp— T()  as  distributed  as  the  age  of  the 
distribution  of  the  time  between  observations,  as  long  as  the  interobservation  times  are  iid.  We  have 
examined  the  simple  case  where  the  interobservation  times  are  constant  at  St;  thus,  tp—  T,  is 
distributed  as  uniform  ( 0,<5 1 ) . 

We  see  that  p^1(  is  calculated  as: 

pfc(<  =  P[  Yr  >  Td  —  tp  and  Ys  >  Td  —  T,  |  on  acquisition  list  at  tp  and  stationary  at  T,  ] 

•  P[  unit  on  the  acquisition  list  at  tp  ]  •  P[  unit  stationary  at  T,  ]  , 


P[  Yr  >  Td  —  tp  |  Td— tp  —  t  ]  P[  Td—  tp  —  t  ]  dt 


E[  Tr  ] 

E[  Tfl  ]  +  E[  Tr  ] 


■f\  \°*  P[Y»>(Td— tp)+(tp— Tj) 
J  1=0  J  i=o 


~ —  *>;  (fp- Tj) — y]  P[Td  —  tp  —  t]  P[tp— T,— y]  dt  dy 


E[S] 

E[  S  ]  +  E[  M  ]  ’ 


[11] 


Example  Case  b 

Suppose  that  the  probability  of  detecting  the  target  during  any  glimpse  period  of  length  6t , 
given  that  it  was  detected  during  the  previous  glimpse  (  pr  ),  is  0.9.  Thus  the  probability  that  we 
fail  to  detect  the  unit  on  the  kth  glimpse  (and  therefore  drop  it  from  our  acquisition  list)  is  equal  to 

(  1  -  Pr  )  Pr*-',  k  =  1,2 . 
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As  a  result,  the  distribution  of  the  retention  time  Tr  is: 

P(  Tr  >  t  ]  =  P(  Kr  >  j-  ],  t  >  0 

where  Kr  ~  geometric  (  1—  pr  ),  and  E[  Tr  ]  =  =  /  — °o~9  =  ^  hours.  Recalling  from  the 

previous  example  that  E[  Ta  ]  =5  hours,  p aCq  =  rr/j?  -  =  0.667. 

1U  "T  0 

Calculating  the  residual  life  (or  age)  of  the  retention  time  is  straightforward.  We  define  an 
integer  m> 0  such  that  m6t  <  y  <  {m+  1)6,  for  any  v  >  0. 

P[  Yr  >  y  ]  =  1  -  -  A-Pr  [  P[  Tr  >  u  ]  du 


1  — Pr  fV1  f fm  +  ijS, 


So  J-J,*  PtT,>u]du  +jtt=m,(P[Tr>u 


1  , r»  771  —  t  , 

=  1 - St  £  Pr  +  (y— m<5()  prm 

°t  k  =  0 

=  1  -  ( 1  —  Pr )  Pp--^Y  -  (y  —  m<5()  prm 

=  1  —  1  Pr  (y  —  mSt)  p rm  ,  m6,  <  y  <  (m  4-1)6,,  m> 0. 

°t  J 

In  order  to  calculate  pa„a,;.  we  need  to  determine  D  =  min{R,  L},  where  R  =  max{  Tr  }  oo 
and  L  =  max{  Td—  tp  }  =  6t  =  I  hour.  Then 

Patat/  =  [  P(  Ar  >  Td  —  tp  |  Td  —  tp  =  t  ]  P[  Td— tp  =  t  ]  dt  •  pac?  •  p3(ay 

_  J  t=o  J 


Pa„a,/  =  0-95  •  Pac,  •  PJ(a„  =  0.95  ■  0.667  •  0.192  =  0.122  . 


We  define  a  time  T,  as  the  last  time  that  the  target  unit  was  observed  prior  to  the  time  of  end 
of  the  planning  process,  tp;  thus  T,  <  tp  <  Td.  Since  6,  is  constant  at  1  hour,  (  tP  —  T,  )  ~ 
uniform!  0,  60  min.  ).  D*  =  min{  L,  Stt  }  =  min{  60  min.,  40  min.  }  =  40  min. 

Ptngaqt  =  Pa„al7  ' 

rso  no 

P[A»>(Td—  tp)  +  (tp—  T;)  |  Td  — tp=t;  tp  — T;=y]  P[Td  — tp  =  t]  P[tp— T,=y]  dt  dy 

J  y=o  J  t=o 

After  some  messy  algebra,  pengage  =  0.122  •  0.039  =  0.005  . 
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Case  c.  The  Target  is  Observed  Continuously  after  Acquisition  wtth  Preplanned  Fire 

This  case  represents  the  situation  when  at  least  one  sensor  can  observe  or  “track”  the  target 
continuously  (or  nearly  continuously)  after  acquisition.  This  may  occur  when  a  potential  target  is 
observed  by  a  fixed  sensor  that  holds  the  target  unit  in  view  (e.g.,  radar).  In  this  case,  we  assume 
that  the  target  remains  on  the  acquisition  list  unit  until  the  sensor  loses  the  target  due  to  movement 
out  of  range  or  into  concealment,  or  until  something  breaks  the  tracking,  such  as  obscuration,  etc. 
The  assumption  of  independence  between  the  acquisition  and  movement  processes  breaks  down  when 
continuous  observation  is  possible.  However,  the  simulation  model  (NUFAM  III)  and  the  acquisition 
data  currently  available  make  this  assumption,  so  we  carry  it  forward  into  our  probability  model 
and  interpret  the  results  accordingly.  CAA  is  in  the  process  of  acquiring  new  data  that  will  enable  us 
to  estimate  the  joint  distributions  directly  (see  section  on  future  research). 


In  this  case,  the  acquisition  is  useful  if  the  target  unit  was  stationary  or  moving  at  the  time  of 
acquisition.  The  probability  that  the  unit  is  on  the  acquisition  list  at  the  planning  time  tp  is: 


P[  unit  on  the  acquisition  list  at  tp  ]  =  pacq 


E[  Tr  ] 

E[  Ta  ]  +  E[  Tr  ] 


[12] 


The  probability  that  the  unit  is  available  for  fire  is: 


P[  YP  >  Td  —  tp  |  on  acquisition  list  at  tp  ]  •  P[  unit  on  the  acquisition  list  at  tp  ],  thus 


Pavai/ 


ro 


t=0 


P[  Y r  >  Trf— tp  |  Tj  —  t p  —  t  ]  P[  Td  — tp  —  t  ]  dt 


ETTri 


E[  Ta  ]  +  E[  Tr  ]• 


[13] 


Since  we  can  track  the  target  unit,  it  need  not  be  stationary  at  the  time  of  acquisition,  but  it 
must  be  stationary  at  the  time  our  fire  planning  is  complete  (tp)  and  remain  at  that  location  until 
detonation  (Td).  The  time  of  detonation  Td  continues  to  be  fixed  during  the  time  of  fire  planning. 
The  probability  of  hit,  pAl(,  is  the  joint  probability  that  the  unit  is  available  for  fire  and  will  be 
stationary  from  the  time  tp  to  Td.  The  times  tp  and  T^  are  preplanned  (independent,  by 
assumption,  of  the  stay/move  process),  so  we  can  simply  multiply  by  p3tay  to  get  the  probability 
that  the  unit  is  stationary  at  time  tp. 


=  P[  Yr  >  Td  —  tp  |  on  acquisition  list  at  tp  ]  •  P[  unit  on  the  acquisition  list  at  tp  ] 

•  P[  Ys  >  Td— tp|  stationary  at  tp  ]  •  P[  unit  stationary  at  tp  ].  [14] 
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P“  =[L  P['’,>T.-1,|TJ-.p  =  llP[TJ-l,  =  tldl].  Err-^rE’[Trl 

'  [L  pnf->T--'.iv>.  =  >iHTi-.,  =  .]d.].  tsrrErri  1151 

Example  Case  c 

Currently,  we  model  a  nearly  continuous  observation  using  the  same  glimpse  model  used  in  the 
previous  sections,  with  a  short  glimpse  period  6,  =  5  min.  If  the  probabilities  of  detection  pa  and 
retention  (1-pr)  are  the  same  as  given  in  the  previous  examples. 

E[  Ta  ]  =  =  25  min.,  E(  Tr  ]  =  ^  =  50  min.,  and  pa<:?  -  0.667  . 

p„.„  =[};:„  p(Y,  >Trtiv.,=,]p[v.,=«]j.]'  • 

r  11  .  1  f  fn  12-ll  " 

=  s,  Pr*  ■  p acg  —  6,  — - j-  ■  pacq  =  0.600  ■  0.667  =  0.100  . 

*=o  J  L  L  Pr_  1  J  J 

no 

Pengage  —  Pavai l  '  P[  V  3  >  Td  —  tp  |  Td~tp  ~  t  ]  P[  Tj  —  tp  =  t  ]  dt 

^  t~°  Ef  S  1 

•  g[  g-j  -;r)  =  0.400  •  0.233  •  0.192  =  0.018  . 

Case  i.  The  Target  is  Observed  Continuously  after  Acquisition  without  Preplanned  Fire 

This  case  represents  the  same  situation  as  case  c  where  at  least  one  sensor  can  observe  or 
"‘track”  the  target  continuously  (or  nearly  continuously)  after  acquisition.  However,  it  is  reasonable 
to  suppose  that  if  we  can  track  the  target  unit  continuously,  we  would  give  the  order  to  fire  when 
unit  stops  for  the  first  time  within  the  fire  execution  period.  Until  now.  we  have  given  the  order  to 
fire  (determined  the  time  Td)  at  time  tp.  In  this  case,  we  define  a  (random)  time  during  the  fire 
execution  period  of  length  L  which  depends  upon  the  time  that  the  unit  stops.  If  the  weapon  can  be 
fired  immediately  upon  detection  that  the  unit  has  stopped,  then  (ignoring  the  time  of  flight  of  the 
round)  Td  =  T j.  However,  it  is  more  reasonable  to  assume  that  some  time  will  have  to  be  spent 
getting  the  order  to  the  delivery  unit,  adjusting  the  aimpoint  to  the  location  at  which  the  unit 
stopped,  confirming  that  the  location  meets  the  criteria  for  employment  of  the  weapon,  etc.  We 
denote  this  time  (Td— T^)  as  77,  and  for  simplicity,  we  assume  that  77  is  constant. 
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If  the  unit  is  stationary  at  time  tP,  Ty=tP:  otherwise.  Ty  equals  the  first  time  the  unit  stops 
between  tP  and  tp  +  (L  — q).  If  the  move  time  distribution  is  not  bounded  from  above  at  some  value 
less  than  the  fire  execution  period  duration,  there  is  a  positive  probability  that  the  unit  will  be 
moving  during  the  entire  fire  execution  period.  The  hit  probabilities  can  be  adjusted  to  account  for 
this. 

Because  the  stay/move  process  is  assumed  independent  of  the  acquisition  process.  T y  is  random 
with  respect  to  the  acquisition  process.  Therefore,  the  probability  that  the  unit  is  on  the  acquisition 
list  at  the  time  Ty  is  simply  pJCl). 

The  probability  that  the  unit  is  available  for  fire  is  still  the  probability  that  it  was  on  the 
acquisition  list  at  Ty  and  retained  between  Ty  and  Td.  However,  if  we  issue  the  fire  order  at  time 
T ,.  we  will  retain  the  unit  between  T y  and  T d  by  definition.  Thus,  the  conditional  probability 
P[  unit  retained  until  Td  |  unit  on  the  acquisition  list  at  Ty  ]  =  1.  As  a  result. 

P Avail  =  *  '  P«<t1 

The  probability  that  the  unit  will  be  stationary  at  the  time  tP  is  simply  pJ(ay.  The  probability 
that  the  unit  stops  within  the  period  (  tP  to  tP  +  (L  —  q)  ).  given  that  it  was  moving  at  time  tP.  is 
equal  to  the  probability  that  residual  life  of  the  move  time.  Ym.  is  less  than  (L  —  q). 

The  probability  of  hit.  p^i(,  is  therefore  the  joint  probability  that  the  unit  is  on  the  acquisition 
list  at  the  stopping  time  Ty.  it  will  be  retained  as  a  target  at  least  until  time  Td.  and  it  will  be 
stationary  from  time  Ty  to  time  Td.  If  Ty  =  tP.  the  unit  was  somewhere  in  its  stay  period  at  time 
tP,  and  the  probability  that  it  will  remain  in  place  until  Td  is  P[  Y,  >  q  ],  uhere  Y.,  is  the  residual 
life  of  the  stay  time.  If  t p  <  Ty  <  L  —  q.  then  the  unit  stopped  during  the  fire  execution  period  and 
Ty  is  the  beginning  of  the  stay  time.  Thus  P[  stationary  at  Tj  ]  tP  <  Ty<  L  — q  J  =  P[  S  >  q  ]. 

p^y  =  P[  Yr  >  Ty  —  Ty  |  on  acquisition  list  at  Ty  ]  •  P[  unit  on  the  acquisition  list  at  Ty  ] 

•  |  P(  Y,>q  |  stationary  at  tp  ]  •  pJ(ay 

+  P[  S  >  q  |  stopped  at  Ty  ]  •  P[  Ym  <  L—  ]  |  moving  at  tp  ]  •  p movc  j 

=  P acq  ■  r  p[  Y ,>T]  }  ■  p)iav  +  P[  S  >  q  ]  •  P[  Ym  L  —  77  ]  •  p move 


L  ^  E[Tr]  __  E[S] 

where  psc,  -  E[  Tfl  ]  +  E[  Tj.  ]  and  P *i«f  ~  E[  S  ]  +  E[  M  ]  ’ 
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Example  Case  d 


From  the  above.  parai/  =  1  •  pac9  =  1  •  0.667  =  0.667 

Suppose  that  the  time  (T^  —  TP).  r /,  is  equal  to  15  min.  The  distribution  of  the  residual  move 
time.  Yin,  is  computed  the  same  way  as  the  distribution  of  Y5. 


P[  Ym  >  y  ]  =  1- 


E[M]  J 


P[  M  >  u  ]  du 


For  0  <  v  <  105  min.,  P[  Ym  >  v  1  =  1  —  . 

■  1  *  ‘  lOo 


P engage  =  Pad,  '  Q  P[  Y,  >  1]  j  •  pj(a!(  +  P[  S  >  T]  ]  •  P[  Yrn  <  L  -  T)  }  ■  pm„( 
=  0.667  •  ^  0.417  •  0.192  +  •  0.808  j  =  0.246  . 


Modeling  Subsequent  Laydowns 

Fp  to  now,  all  of  the  development  has  been  centered  upon  a  single  nuclear  laydown  per  side, 
which  involves  a  single  planning  period  and  single  reference  point  tP  at  which  we  determine  the  unit 
acquisition  and  movement  states  and  the  associated  residual  life  distributions.  If  we  wish  to  model 
subsequent  laydowns,  we  will  need  to  make  determinations  of  the  unit  states  at  some  later  planning 
time,  say  t  P* .  where  tP*  >  tP.  If  tP*  and  tP  are  sufficiently  far  apart  (this  is  a  matter  of  judgement 
—  generally  at  least  the  expected  length  of  several  transition  cycles),  then  the  state  determination  for 
each  unit  at  time  tP*  can  be  treated  as  independent  of  the  state  determination  at  time  tp.  If  tP*  is 
close  in  time  to  tP,  the  state  determinations  will  be  dependent  and  this  dependency  should  be 
represented  in  the  model.  For  a  discussion  of  these  points,  along  with  a  suggested  approximation  for 
this  dependency  relationship,  see  Youngren  [1989]. 

Applications  of  this  Methodology  at  CAA 

The  US  Army  Concepts  Analysis  Agency  uses  the  methods  described  previously  to  determine  (he 
probability  that  relocatable  target  units  may  be  acquired  and  hit  by  tactical  nuclear  weapons  at 
times  when  running  the  detailed  simulation  model  NUFAM  III  is  not  possible  or  desirable.  These 
methods  permit  an  estimation  of  the  impact  of  varying  input  data  (for  example,  different  values  for 
the  average  stay  and  move  times  for  units)  on  NUFAM  output,  without  running  the  NUFAM 
model.  These  same  methods  can  be  incorporated  into  NUFAM  III  to  replace  the  explicit  simulation 
of  the  acquisition  and  movement  events.  Such  a  replacement  eliminates  over  2  million  simulation 
events  per  NUFAM  excursion  without  a  loss  of  relative  accuracv  in  the  results  of  interest. 

1 
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Future  Research 

To  eliminate  the  necessity  of  assuming  independence  between  acquisition  and  movement,  and  to 
improve  our  estimates  of  acquisition  and  the  probability  of  hitting  the  target.  CAA  has  an  ongoing 
data  collection  effort  with  the  US  Army  Systems  Analysis  Agency  to  try  to  determine  the 
probabilities  of  acquisition  and  retention  under  multiple  factors  such  as  stationary/moving;  posture; 
weather;  obscuration;  sensor  system  deployment,  etc.  Detailed  sensor  models  will  be  used  to  estimate 
the  distribution  functions  for  Ta  and  Tr  conditioned  on  these  factors.  The  Agency's  target 
acquisition  model  will  be  modified  to  vary  the  factors  according  to  their  underlying  distributions  in 
order  to  estimate  the  unconditional  joint  distributions  of  Ta,  Tr.  S.  and  M  for  the  scenario  of 
interest.  These  distributions  will  be  used  to  improve  our  estimates  of  acquisition  and  engagement 
and  to  refine  the  probability  model  to  account  for  the  dependencies  between  movement  and 
acquisition. 

The  probability  models  discussed  in  this  paper  will  eventually  be  incorporated  into  a  new 
stochastic  model  of  a  theater-level  nuclear  exchange  called  NEMESIS,  which  will  replace  NUFAM 
III.  NEMESIS  will  combine  the  probabilistic  description  of  our  uncertainty  about  the  high  resolution 
processes  of  target  acquisition  and  movement  given  in  this  paper  with  a  probabilistic  description  of 
our  uncertainty  about  the  locations  of  the  potential  target  units,  providing  increased  resolution  in 
space  as  well  as  time  over  theater  models  such  as  FORCEM.  The  NEMESIS  model  may  either  be 
used  as  a  stand-alone  model  of  theater  nuclear  warfare  (in  a  manner  similar  to  NUFAM  III)  or  as  a 
means  of  generating  realizations  of  possible  nuclear  exchange  outcomes  for  input  into  the  FORCEM 
model.  The  latter  use  permits  us  to  model  the  effect  of  different  possible  nuclear  exchange  outcomes 
on  subsequent  conventional  battle  in  FORCEM. 
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